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Abstract 

In this paper, we consider the Stokes equations and we are concerned with the inverse problem of 
identifying a Robin coefficient on some non accessible part of the boundary from available data on the 
other part of the boundary. We first study the identifiability of the Robin coefficient and then we establish 
a stability estimate of logarithm type thanks to a Carleman inequality due to Bukhgeim |^ . 
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1 Introduction 

Let us consider some open Lipschitz bounded connected set domain of R'^, d > 2. We assume that the 
boundary F = dfl is composed of two parts Fq and Fg such that Fe U Fq = F and Fg n Fo = (Figure [I] gives 
an example of such a geometry) . We introduce the following boundary problem: 

utit,x) - Au{t,x) + Vp{t,x) = 0, yxen,yt>o, 
y-u{t,x) = 0, yxefi,yt>o, 

Vu{t,x)-n{x)~p{t,x)n{x) = git,x), Va;GFe,Vt>0, (1.1) 

Vuitjx) ■ n{x) ~ p{t,x)n{x) + q{x)u{t,x) — 0, Va;eFo,Vi>0, 
u{0,x) — uq{x), \fx Q fl. 

Notice that we assume that the Robin coefficient q defined on Fg only depends on the space variable. Our 
objective is to determine the coefficient q from the values of u and p on Fg. 

Such kinds of systems naturally appear in the modeling of biological problems like, for example, 
blood flow in the cardiovascular system (see 201 and |23]) or airflow in the lungs (see (Sj). The part of 
the boundary Fg represents a physical boundary on which measurements are available and Fq represents an 
artificial boundary on which Robin boundary conditions or mixed boundary conditions involving the fluid 
stress tensor and its flux at the outlet are prescribed. 

Similar inverse problems have been widely studied for the Laplace equation [1], [10], [H], [H] 
and |22| . This kind of problems arises in general in corrosion detection which consists in determining a Robin 
coefficient on the inaccessible portion of the boundary thanks to electrostatic measurements performed on 
the accessible boundary. Most of these papers prove a logarithmic stability estimate ( [2], [1], [IQ] and [12j). 
S. Chaabane and M. Jaoua obtained in JJJ both local and monotone global Lipschitz stability for regular 
Robin coefficient and under the assumption that the flux g is non negative. Relaxing this constraint, they 
obtained in a logarithmic stability estimate. More recently, E. Sincich has obtained in [2^ a Lipschitz 
stability estimate under the further a priori assumption that the Robin coefficient is piecewise constant. To 
prove the stability estimates, different approaches are developped in these papers. A first approach consists 
in using the harmonic functions properties (see [2], [10]). A characteristic of this method is that it is only 
valid in dimension 2. Another classical approach is based on Carleman estimates (see [1] and [H]). In [3], 
the authors use a result proved by K.D. Phung in [19J to obtain a logarithmic stability estimate which is 
valid in any dimension for an open set fl of class C°° . This result has been generalized in [S] and to C^'^ 
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and Lipschitz domains. Moreover, in 0], the authors use semigroup theory to obtain a stabiHty estimate in 
long time for the heat equation from the stabiHty estimate for the Laplace equation. 



In this article, we prove an identifiability result and a logarithmic stability estimate for the Stokes equations 



with Robin boundary conditions (1.1 1. The paper is organized as follows. 

The second section contains preliminary results on the regularity of the solution. For this purpose, 
we introduce Besov spaces: they appear as the natural space to which the Robin coefficient q belongs. 

In the third section, we are interested in the identifiability of the Robin coefficient q. Under some 
regularity hypotheses and using the theorem of unique continuation for the Stokes equations proved in |15j . 
we prove that if two measurements of the velocity are equal on the boundary (0,r) x Fg, then the two 
corresponding Robin coefficients are also equal. 

Section [3] corresponds to the main part of our article. The results of this section are only valid 
in dimension 2. We prove a logarithmic stability estimate, first for the stationary problem and then for 
the evolution problem. To do this, we use a Carleman inequality due to Bukhgeim which is only valid in 
dimension 2 (see [S]). The stability estimate for the unsteady problem is deduced from the stability estimate 
for the stationary problem thanks to the semigroup theory. We end Section [4] by concluding remarks and 
perspectives to this work. 

When we are not more specific, C is a generic constant, whose value may change and which only 
depends on the geometry of the open set and of the boundaries Fg and Fq. 

We are going to start with some preliminary results which will be useful in the subsequent sections. 

2 Preliminary results 
2.1 Besov spaces 

We introduce the following functional spaces: 

V ^{v e H^{nf/V ■v = Oonn}, 

and 
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We also need to introduce the space to which the Robin coefficient q belongs: we are going to assume enough 
regularity on q in order to give a sense in suitable spaces to the trace product qu. We refer to |16| for more 
general spaces. Let seR, l<p<oo and: 

where S'iW^) is the space of tempered distribution on and J-'w is the Fourier transform of w. Equipped 
with the norm 

iMs^.^iR") 11(1 + icn*-7^wii^p(Kd), 

i?s,p(IR'*) is a Banach space. We observe that Bs^2{^'^) = H''{R'^). In the following, we will need the following 
properties: 

Proposition 2.1. Let s > 0. We have Bs.i{k'') ^ C°{R''). 

Proof of Proposition\2.1\ Let u G i3s,i(IR'*). We start by proving that J^u e L^(IR'*): 
and from the Cauchy-Schwarz inequality, we have: 
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iJ"«(OMe< ii(i + ien*-^"iiLMK-) 



<Ch|lB,,(R^). (2.1) 



Since for almost every a; G R'^, 



we deduce the continuity of u from the Lebesgue continuity theorem and according to (2.1 1, we have: 

ll"llco(R<i) < C'll'"||B,_i(Rd). 

□ 

Using a local map and partition of unity, we build i3s,p(r) from Bg.piR"^^^) in the same way that 
H'^{T) is built from H''{R'^^^) (see ^7\). From the above proposition, we deduce the following corollary: 

Corollary 2.2. Let s > 0. We have i?s^i(r) ^ C'^(r) and the injection is continuous. 

Proposition 2.3. Let s e R+ , u E iJ^(r) and v e Bs,i{r). Then uv G H''{T) and 

\\uv\\hs(t) < (2^)^-1 H"Hg°(nlMlB.,i(r)- 
We refer to [T3] for a proof of this proposition. 

2.2 Regularity of the stationary problem 

First, we are interested in the stationary case: 



-Au + Vp = /, in n, 

V • u =0, in f], 

Vu-n — pn ~ g, on Fg, 

Vu-n — pn + qu — 0, on Fg. 



(2.2) 



Let g e H-2{T^y and v G H2(TeY\ we denote < g,v >_ i i the image of v by the linear form g. 
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Proposition 2.4. Let a > 0, f e L'^{^Y, g € H-'siT^)'' and q e L°°(To) such that q > a ohTq. System 
(2.2 I admits a unique solution {u,p) G x Moreover, there exists a constant C{a) > such that 



IMmin)" < C'(")(ll.9ll^-i(p_^)<i + \\.f\\L-Hny) 



(2.3) 



Proof of Proposition 2.4 The variational formulation of the problem is: find u £ V such that for every 

veV, 



Vu : Vu + / qu ■ V =< g,lY^v >_i 1 Y + f-v 



We note V(u, v) e 



and Vf £ V, 



aq{u,v) — / Vu : Vw + / i 



qu ■ V, 



(2.4) 



We easily verify that Oq is a continuous symmetric bilinear form. Since q > a > 0, according to the 
generalized Poincare inequality, the bilinear form Oq is coercive. On the other hand, Li is a continuous 
linear form on V. Thus we prove the existence and uniqueness of m G ^ solution of equations (2.2) using the 
Lax-Milgram Theorem. We obtain simultaneously estimate (2.3 1. We prove the existence and uniqueness of 
p g L^(f2) in a classical way, using De Rham Theorem (we refer to [7] for the case of Neumann boundary 
condition) . □ 

Let us recall existence and regularity results for the Stokes problem with Neumann boundary 
condition proved in 7J. 

Proposition 2.5. Let /c e N. Assume that is a bounded and connected open set in R'* of class C'^^^'^. 
We assume that: 



Then the solution {u,p) of the problem: 



-Alt + Vp 

V • u 
Vu • n — pn 



f in Q, 
in J7, 
h on r. 



belongs to H'^^^[il)'^ x H^^^{Vl) and there exists a constant C > such that: 

WAm+HU)'^ + Ilp|lff^-+i(i2) < C(||/i||^i+fc^^^^ + \\f\\H''(n)i)- 

From the previous proposition we deduce the following result: 

Proposition 2.6. Let A: G N. Assume that is a bounded and connected open set in R'^ of class C'^^^'^. 
Let a > 0, M > 0, / e H''{n)'', g e Hi+''{T^Y and q e Bi+j. ^(ro) such that q > a on Tq. Then the 

solution {u,p) of system ( |2.2[ ) belongs to i/'^+^(il)'^ x ij'^+-'^(il). Moreover, if\\q\\Bi^i^ ^{To) — -^^j there exists 

a constant C{a,M) > such that 

MH^+^ny + MH^+^n) < ^(a, Af)(||5||^,+ i^^^^_^ + UWh^u)")- 

Proof of P roposition \2. 6\ Let us prove the result for fc = 0. Let h = — qulrn + g'^r^- According to Propo- 
sition 2.4 u belongs to H^{n)'^. Thus, we obtain from Proposition 2.3 that qu € H^iT^)''- , which implies, 

with fc = we obtain that 



2.5 



since g e Hi{TeY and Pg n Fo =0, that h e iJ3(r)''. Using Proposition 
{u,p) e H^iny^ X H^{n) and: 
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But, we have from Proposition |2.4[ that: 



and since by hypothesis, H^Usi (Fq) < we obtain: 



Thus we obtain the result for fc = using the inequahty of Proposition |2.4[ We then proceed by induction 
to prove the result for any /c G N. We emphasize that the key argument is tha t if m G then 
u G _ff'^+2 (Fo)'' which implies that g u^-p^ £ _ff 2+'''(ro)'^ thanks to Proposition 
regularity result given by Proposition 2.5 and conclude. 



2.3 



Thus, we can apply the 

□ 



2.3 Regularity of the evolution problem. 



Concerning the initial problem (1.1), we can prove, using Galerkin method, the following regularity results. 
For completeness, the proof of Theorem |2.7| is given in the appendix. 

Theorem 2.7. Assume that fl is a bounded and connected open set in R'^ of class C^-^ Let T > 0, M > 0, 

a > and uq G V. We assume that g G H^{0,T\ (T^)'^)^ q g Bi i{To) such that \\q\\Bi (To) 1^ ^'^ ^'^'^ 

2 ^ 2 

q>a on Tq. Then problem ([Tl]) admits a unique solution {u,p) G L'^{0,T; H'^inY) f] H^{0,T; L'^ifl^) f] 
i°°(0,T;y) X L^(Q,T-H^{n)). 

The following corollary will be useful when we will prove logarithmic stability estimate for the 
evolution problem (1.1). 

Corollary 2.8. Assume that fl is a hounded and connected open set in R'^ of class C'^'^ . Let M > 0, 
T > 0, a > and Uq e H^{ny f] H. We assume that g G iJ2(0, T; i (Fg)''), q G B3 ^^{Tq) such that 
W^Wbs j(ro) ^ q > a onTo. Then, problem ( 1.1 ) admits a unique solution {u,p) G L°°(0, T; H^{il)'^)r\ 

H\0,T; H'^in)'^) n H^{Q, T; L'^{n)'^) x L°°(0, T; H^{n)) D H\Q, T; H\n)). 



The proof of the previous corollary consists in applying Theorem 2.7 to ut- Let us prove it. 
Proof of Corollary\2.8\ Let {u,p) the solution of (1.1). Let us consider the following system: 



vt — Av + Vt 


= 0, 


in (0, T) X n, 


V • V 


= 0, 


in (0,T) X n, 


Vu ■ n — Tn 


= 9t, 


on (0,T) X Fe, 


Vu ■ n ~ Tn + qv 


= 0, 


on (0,T) X Fo, 


v{0) 


= Amo-Vp(O), 


in fl, 



(2.5) 



where p{0) G H'^{VlY is defined as the solution of the following elliptic boundary problem: 

Ap(0) = 0, in f7, 

p(0) = {Vuq- n) ■ n- g{G) ■ n, on Fe, 
p(0) = (Vwo • n) • n + g(0)Mo • n, on Fq. 



2.6 



According to Theorem Et) we obtain that (w,t) belongs to L'^{Q,T] H'^iVlY) n H^{0,T; L'^{VLf) n 
L°°(0,T;y) X L^{0,T;H^n)). Remark that {ut ,Pt) is solution of equations ( |2.5| ) in the distribution sense 
on (0,r). Thus, by uniqueness, (w,t) = {ut,pt)- Then, since q G -Bi liTe), we deduce from Proposition 
that the linear map: 

K(t),.g(i)) ^ {u{t),pit)) 

is continuous. Since {ut,g) G L°°{0,T;V) x L°°{0,T; H^^e)'^), we deduce that (u,p) G L°°{0,T; H^inf) x 
L°°{0,T;H^{n)). □ 
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3 Identifiability 

3.1 Unique continuation 

We start by recalling a unique continuation result for the Stokes equations proved in [15] . 

Theorem 3.1. Let be an open connected set in R^, d> 1. We note Q = (0,T) x and let O be an open 
set in Q . The horizontal component of O is 

C{0) = {{t, x) e g, 3x0 e (t, xo) e O}. 

Let {u,p) e L^{0,T;Hl^^{n) f x Lf^^{Q) be a weak solution of 

Ut-Au + Vp = 0, in{0,T)xn, 
V-u = 0, m(0,r)xf], 

satisfying u = in O then u — and p is constant in C{0). 

We easily deduce the following result from the previous result. It will be very useful in the next 
subsection. 

Corollary 3.2. Let 6 > , xq Cz T and r > such that j = {to — 6, to + S) x {B{xo, r) D T) is an open set in 
(0,r) X r. Let {u,p) e L^{Q,T-H^{nY) X L^{0,T;H^{n)) be solution of: 

ut- Au + \/p 0, in (0, T) x n, 
V-u = 0, inlo,T)xn, 

satisfying u = and Vm • n — pn = on 7. Then u ~ and p ~ in (to — d, to + 5) x fi. 
Proof of Corollary \3.S\ We extend u and p by on (^o — 5,to + 6) x {B{xo, r) n fi^): 

^ P " \ in {to ~S,to + 5)x {B{xo, r) n n^) 

and we denote VL ^ n\J B{xo,r). Let us verify that {u,p) e L'^{0,T; H^{n)'^) x L'^{0,T; L'^{n)) is stiU a 
solution of the Stokes equations in il. Let v e T>{n)'^. We check by integration by parts in space that almost 
everywhere in i € (io ~ S,to + S): 

ut ■ V + / Vu • Vw — / pV ■ V — 0. 
n Jn Jn 



Moreover V ■ u = in {to — 5,to + S) x fl. Therefore, we can apply Theorem 3.1 to {u,p): {u,p) = (0, 0) in 
{to — S,to + S) X which implies that u — and p is constant in {to — 5,to + 5) x il. At last, the fact that 
Vm • n — pm = over 7 implies that p — va {to — to + 5) x Q. □ 



3.2 Application 

Proposition 3.3. Let T > 0, a > 0, x^ eT^, r > 0, g e H'^{0,T; {TeY) non identically zero, uo €z V 



and Qj G Bi i{To) such that qj > a onTo for j — 1, 2. Let {uj,pj) be the weak solutions of (1.1 ) with q — qj 
for j — 1, 2. We assume that ui = U2 on (0, T) x {B{xe, r) f] Fe). Then qi = q2 on To- 



Proof of Proposition \3.3\ We are going to prove Proposition |3.3| by contradiction: we assume that qi is not 
identically equal to 52- 

We have {uj,pj) e L^{Q,T- H'^{VLf) x L^{0,T: H^{n)) for j = 1,2 thanks to Theorem 



2.7 



We define 



u = ui — U2 and p — pi ~ p2- Let us notice that is the solution of the following problem: 

ut~ Au + Vp = 0, in (0, T) x fi, 

V-u =0, in(0,T)xr2, 

Vu-n-pn =0, on(0,T)xre, 

Vu ■ n — pn + qiUi — (72^2 = 0, on (0, T) x Fq. 
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By hypothesis, u = and Vu • n — = on (0, T) x {B{xe,r) H Fg). Thus, according to Corollary 3.2 
ui — U2 and pi — p2 in (0, T) x 57. Consequently, we deduce from 

Vui ■ n — Pin + qiUi ~ 0, on(0,r)xro, 
Vui ■ n — Pin + q2Ui = 0, on(0,T)xro, 



that 



Mi(gi-g2)-0on (0,r)xro. 



(3.1) 



By hypothesis, qi is not identically equal to q2- Thanks to Corollary 2.2 qi and q2 are continuous. Then, 
we can find an open set k d Tq with a positive measure such that: 

(<?! - q2){x) ^ 0, Va: € K. 



Equation (3.1) implies that ui = on (0,r) x k and then ui is the solution of 



uit-Aui+Vpi = 0, in(0,r)xrj, 

V-ui 0, in(0,r)xf7, 

ui = 0, on(0,r)xK, 

Vui • n — Pin = 0, on (0, T) x k. 



Applying again Corollary 3.2 we obtain that ui ~ Q and pi = in (0, T) x f2. This is in contradiction with 
the assumption that g is non identically zero. □ 



4 Logarithmic stability estimates 



In this section, we assume that d — 2 and that the open set £7 C IR is of class C ■ in order to obtain regular 



solutions of problem (2.2) from Proposition 2.6 We are going to prove logarithmic stability estimates using 



a Carleman inequality which is stated in Lemma [4~T| First, in Theorem |4.3[ we state a logarithmic stability 
estimate for the stationary problem. Then we deduce from this Theorem two logarithmic stability estimates 



for the evolution problem using the analytic semigroup theory. These estimates are given in Theorem 4.16 



when the flux g is stationary and in Theorem 4.18 when g depends on time 



4.1 Carleman inequality 

Let us first state a Carleman inequality proved by Bukhgeim [9]: 
Lemma 4.1. Let e C'^ (Tt) . We have: 

f (A*|m|2 + (A*- l)|Vu|2)e* < / |Au|2e*+ / V^iJ ■ nHuf + \Vuf + 2\dr\Vuf\)e^ (4.1) 
Jo Jo Jr 

for allue C^{n). 

The proof of this result, which is only valid in dimension 2, uses computational properties of 
function defined on C (in particular, the fact that 4:dzdz — A). 

Remark 4.2. The result is still true for u £ H^{fl). Indeed, for all u G H'^{VL), there exists (M„)„g[M G 
C2(n)'^ such that 

Un^uinH^{n). (4.2) 
We can apply Lemma\4.1\ to u„, for all n £ N. Let us prove that: 



lim / V*-n|9^|Vw„ne* = / V* • nia^lVwHe*. 

n^oo J J. _/p 



(4.3) 



Note first that Jp V^* • n|9r|VMp|e* has a meaning for u G H^{il): 



^ Vvl/ . n\dr\Vuf\e^ < 2||vl/||c,(n)l|e*|lco(o) ' l^^^^l) < ^■ 
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We have: 



IV^- • n\ WdJVi 



2| I „* 



<CmcHn)\\^''\\coin) E / l^^" 



According to (4.2), t/ie sequence {dijUn)n(£\ti converges in L'^(T) towards dijU and \\dijUn\\L'^(^r) *s bounded by 
a constant independent ofn. Ttien, equality (4.3) follows from (4.2). 



(4.4) 



4.2 The stationary case 

For the stationary problem: 

-Aw + Vp =0, in n, 

V - u =0, mQ, 

Vu-n—pn — g, on Fg, 

\7u-n — pn + qu = 0, on Fq. 

we have the following stability estimate. 

Theorem 4.3. Let a > 0, Mi > 0, A'h > 0, {g,qj) G Hi{Te)^ x ^(Fo) for j = 1,2 such that g is not 
identically zero, \\g\\ 5 < Mi, qj > a on Tq and HqjUs^ (Fo) — -^2- We note {uj,pj) the solution of 

(4.4) associated to qj for j = 1,2. Let K be a compact subset of {x € Fq ui 7^ 0} and let m > be a constant 
such that \ui\ > m on K . Then there exist positive constants C{m, Mi, a) and Ci such that 



Iki - q2\\mK) 



< 



C(m, Ml, Af2, a) 



In 



Ci 



■"i-"2|| £,2(r^)2 + ||V(tii-U2)-n|| i2(r^)2 + ||pi-p2|| i,2(r^) + l|V(pi-p2)-n||_c2(r 



Remark 4.4. Note that, thanks to Corollary 3.2, we ensures that {x € Fq ui 7^ 0} is not empty. Then, 
thanks to the continuity of ui, we obtain the existence of a compact K and a constant m as in TheoremU- 



Remark 4.5. In the same kind of inequality is proved for the Laplacian problem with Robin boundary 
conditions under the hypothesis that the measurements are small enough. Here, we free ourselves from this 
smallness assumption on the measurements. 

Remark 4.6. Remark that, thanks to the boundary condition on Fg, inequality in Theorem \4-.3\ can be reduced 
to: 

,1 ^ C(m,Mi,M2,a) 

ll^l - q2\\L^K) < 

'in 



ll'"i-«2|li2(r^)2 + |lpi-P2|| t2(r^) + l|V(pi-p2)-n|li2(r^) , 

Remark 4.7. Comparing this result with Proposition \3. 3t we can emphasize some differences. In Proposi- 
tion's^ we only need that ui — U2 on j where j is a part of the boundary included in Fg in order to recover 
the Robin coefficient everywhere on Tq. However, the proof uses the fact that the constraints are equal on 
7. Here, we need to have measurements on the whole part of the boundary Fg, and the constraint is divided 
into two terms: ||V(ui — U2) ■ ^Ili2(p y in one hand and \\pi "P2|Il2(p ) in the other hand. Moreover, we 



have an additional term: ||V(pi 



P2) 



L2(re)- 



Let us begin by proving this intermediate result which gives us a logarithmic estimate of the traces 
of u, Vii, p, Vp over Fq with respect to the ones over Fg. 

Lemma 4.8. Let {u,p) € iJ^(ri)^ x H^{^) be the solution in fl of 



-Au + Vp 
V • u 



0, 
0. 



We assume that there exists A > such that 

\W\\H'i(ny + \\p\\m{n) < A. 



(4.5) 
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Then there exist C{A) and Ci such that: 
ll"lli^(ro)2 + ||Vu||i2(ro)4 + ||p||i2(ro) + ||Vp||L2(ro)2 < 



In 



C{A) 

Ci 



■"||£,2(r^)2+||Vn-n||_j_2(r^)2 + ||p||£,2(r„i + II^P'" 



I,2(re) 



Proof of Lemma 4-8 The proof is based on the Carleman inequaUty of Lemma 4.1 for an appropriate 



choice of Note that we wiU use (4.11 twice: one time for the velocity u and one time for the pres- 



sure p. The weight function ^ is chosen in order to estimate the traces over Fq with respect to the ones on Fg. 
Step 1: choice of Vt. 

We choose 5* as in [TJ]. There exists ^'o £ C^(fi) non identically zero such that: 

A*o = in fi, = on Fq, fo > on Fe, V^fQ • n < on Fq. 
Indeed, let x ^ C^(F) such that 

X = on Fq, X > on Fe, 
and X noil identically zero on Fg.The boundary value problem : 

f A^-o = 0, in n, 
\ ^-0 = X, onF, 

has a unique solution vj/g g C^(r2). Note that is not constant because x is non identically zero. So, from 
the strong maximum principle, ^'o > in f2. According to Hopf Lemma, we have V^o ■ < on Fg. 

Let A > 0. Denote *i G C^iU) the unique solution of the boundary value problem: 



^-1 



A, in 
0, on F. 



From the comparison principle and the strong maximum principle, we have ^'i < in f2. Moreover, 
according to the Hopf Lemma, we have V^i • n > on F. 

Let us consider 4* = + s^'o for s > 0. To summarize, the function has the following proper- 
ties: 

A* = A in rj, * = on Fq, > on F^, and sV*o ' " < • n < V^-i • n on Fq. 



Step 2: We first apply Lemma 4.1 to u. Using the fact that Au — Vp, we have: 



/ (A*|w|2 + (A^- l)|VMne* < / |Vp|2e*+ / ■ ni\u\^ + \Vu\^ + 2\dr\Vu\^\)e'^ 
jQ Jq Jr 



(4.6) 



We then apply once again Lemma 4.1 to p, it yields: 

/ (A^-IpI^ + (A^ - l)|Vpne* < / lAppe* -f / V* • + [Vpp -t- 2\dr\\/p\^\)e^ . (4.7) 

Jn Jn Jr 



We have Ap = div(AM) = hence |Appe* = 0. We now choose A > 2. By summing up inequalities (4.6) 
and (4.7) and eliminating the integrals over D, in the left hand side which are positive terms, we obtain: 



J V^^-ndwp-f |Vup + 2|9,|Vun)e*+ J V* • 7i(|p|2 -f |Vpp + 2|9,|Vpn)e* > 0. 
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We now specify the dependence with respect to s. We denote 6 — niinp^i [V^'o ' ^Ve note that on Tq? 
e* = 1. Consequently: 



+ 2 [ V* • n(|9,|Vpn + \dr\S7u\^\) + 2 f V* • n(|a,|Vpn + la^jVupDe* 

+ ^ • + |Vu|2 + + iVp^e* > 0. (4.i 

Let us study each of the terms. We have: 

VvPi . n{\u\' + |Vup + bp + iVpp) < C(||u||^3(s,). + \\p\\l.(n))- 



Moreover, since V^* • n < V^*! • n on Fq and using hypothesis (4.51, we obtain 



2 / . ni\dr\Vp\'\ + |9.|Vwh) < C{A)i\\u\\ 



Since, on Fg, jV^* • n\ < sC for s > 1 , using hypothesis (4.5) and Cauchy-Schwarz inequahty, we obtain: 



2 J V^-n{\dr\Vp\'\ + \dr\Vu\'\)e'' <sCiA)l^J (|Vpp + |Vup)e 



Similarly, for s > 1 we have: 



V* • n(bp + \Vu\' + br + |Vpp)e* < Cs / (bP + + \p\ + \^P\ )e 



Note that e* depends on s over Fg. Hence, reassembling these inequalities, inequality (4.8) becomes: 
bp + iVzip + bp + iVpp) < CiA)K, + -i\\u\\Un)^ + \\p\\Un)), 



where Ks^ /p JbP + l^up + bP + |Vpp)e* + (^/p_^ (| Vpp + |Vup)e^* j . Remark that, thanks to classical 
interpolation inequalities (see [1]), there exists C > such that for all / G i?^(Fe): 



Applying the previous inequality and the assumption (4.5), there exists C{A) > such that: 

J^JVu-t\^ <C{A)\\u\\mr.)., and /pjVp • rp < C(A)|blU.(r.). 
In order to precise the dependence with respect to s of if^, we denote: 



6 = 



bP + |VM-np + |Vp- 



(4.9) 



(4.10) 



We obtain, using the fact that Vw = Vu ■ n + Vu • r on Fg and inequality (4.9), that there exists C{A) > 
such that: 

Ks<e'''\ I bp + |Vup + bp + |Vpp + (^y |Vj5p + |Vup) ] <CiA)e''%6 + Si +Si), 



where k = maxp^ ^'o- To summarize, using again assumption (4.5), we have that there exists C{A) > such 
that: 

Ks <C{A)e'''Si. (4.11) 
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Hence we get for all s > 1: 

C 

To 



Remark that this inequality is trivially verified for < s < 1 by continuity of the trace mapping. To 
summarize, we have proved: 



+ \Vu\^ + bp + iVpp) < C{A) [e'^'S'i + i j , Vs > 0. 

We now optimize the upper bound with respect to s. We denote /(s) = e'^^Ji + K Let us study the function 
/ in R:;. We have: 

lim^^o I{s) = +00, 
lims^oo /(s) = +00. 

So since / is continuous on 1R+, f reaches its minimum at a point Sq > 0. In this point, 

/'(so) ^0<^6^ ^ thus /(so) 

Hence: 

' (|up + |v^.p + H2 + |vpn< ^ 

where /3 = 1 if so > 1 and /? = 2 otherwise. But: 



1 


1 


fcSQ 


\ — 

So 






a- 







4 - A^so'e'^^" < fce('=+2)^«, 
(5* 



that is to say: 

1 ^ fc + 2 



so In ( 

\kSi J 

In the same way, when sq < 1 , we obtain: 



1 ^ k + l 



4 I 



( I \ 



Using the fact that Zn ^xa^ — ^ln{x) for all a; > and remembering the definition (4.10) of d, the desired 
result follows. □ 



Let us now prove Theorem |4.3[ 

Proof of Theorem \4-3\ We have on Fo : 

(g2 - gi)ui = g2(ui - U2) + V(ui - U2) ■ n - (pi - p2)n, 

which leads to 

\\qi - q2\\L^(K) < C{m,M2) (||ui - U2||L2(ro)2 + l|V(wi - U2) ■ n||L2(r„)2 + \\pi -p2||L2(r„)) 



Let j = 1,2. According to Proposition 2.6 the solution {uj,pj) of problem (4.4| associated to qj belongs to 
H^{VlY X H^{Vl) and moreover there exists C{a,Mi,M2) > such that 

||wi - U2||i/4(Q)2 < C{a,Mi,M2) and \\pi -p2\\H^n) < C(a, Afi, M2). 

Consequently, we can apply Lemma |4.8| and we obtain: 

C{a,m,Mi,M2) 



hi - q2\\L^K) 



< 



In 



"i-"2|| i,2(r^)2 + ||V(tii-U2)-n|| i2(r^)2 + ||pi-p2|li,2(r^) + l|V(pi-p2)'n||j:,2(r^) 



□ 
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4.3 Evolution problem 

In order to use semigroup properties, we begin by introducing the Stokes operator associated with the Robin 
boundary conditions on Tq. 

4.3.1 Properties of the Stokes operator 



We recall that the bilinear form is defined by (2.4 1. 

Definition 4.9. We define the set 'D{Aq) as follows: 

V{Ag) = {ue V/3C > 0,V« e V, \a,{u,v)\ < CMmny}, 

and the operator Aq : ^{Aq) d H ^ H by: 

Vu e VlAq), aq(u,v) = (^qM, u)i2(n)2,Vu e V. 

Proposition 4.10. Let a > and q G L°°{Tq) such that q > a almost everywhere on Tq. The operator Aq 
has the following properties: 

1. Aq G C{'D{Aq), H) is invertible and its inverse is compact on H. 

2. Aq is selfadjoint. 

As a consequence, Aq admits a family of eigenvalues (pq 

Aq<j)q = \q<j)q with < < < ... < and lim A^ = +00, 

which is complete and orthogonal both in H and V . 

Proof of Proposition \4-l(^ It relies on classical arguments for which we refer to [5j or [21] . □ 

1 

Corollary 4.11. The operator Aq is an isometry from (V, || ||_ffi(o)) lo {H,\\ \\L^(fi)). 

According to the min-max Theorem, since aq(w, u) > aa{u, u) for all u £ V, we have the following 
lower bound: there exists fi > such that for alH > 1 

A^ > (4.12) 

Proposition 4.12. Let a > and q G L°°{Tq) such that q > a almost everywhere on Tq. The operator 
—Aq generates an analytic semigroup on H . This analytic semigroup is explicitly given by: 

e-*^'/ = E^"'''('^U)i^(o)^<' (4.13) 
i>i 

for all f e H. 

Proof of Proposition \4- 1S\ It follows from the construction of the operator Aq. We refer to 18J and 14] for 
details. □ 

Proposition 4.13. Let fc G IM. We assume that q G Bi^f, liTa) is such that on Fq, q> a. Then for each 
f <E H f] H^{^)'^, there exists u G -ff'^+^(0)^ solution of AqU — f if and only if there exists p G H''^^{fl) 
such that {u,p) is solution of the following problem: 

—Au + Vp = / in n 

V-u — in , , -, ,^ 

Y7 n r (4-14 

Vu ■ n ~ pn — U on L e 

Vm ■ n — pn + qu = on Fq 

Moreover, if we assume that there exists M > 0, ||(z||bi^^ ^(Fq) ^ j then there exists C{a,M) > such 
that ||w||^2+fc(o)2 < C(a,M)||/||^fc(n)2. 
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Proof of Proposition \4-.l^ It follows from the construction of the operator Ag and from Proposition 2.6 □ 
Corollary 4.14. Let fc e N* and q e Bi^2{k-i),ii^o) such that q > a ouTq. Then V (A'^j) '-^ H^'' (n)^ n H . 

Proof of Corollary \4-14 For fc = 1, it is clear. Take now fc = 2. Let u e V{Al). We have 

A^v - f^l AqU = V 

Agv = f 

But V g VlAq) C H^{Vl)^ n by assumption, so m g H^{HY H H thanks to the regularity properties of the 
solution of the Stokes problem summarize in Proposition |2.6[ We conclude by induction on fc. □ 



Remark 4.15. Let us remark that'D{Aq) is not equal to H^{fl)^nH : it comes from the boundary conditions. 



4.3.2 The flux g does not depend on t 



In this paragraph, we consider the evolution problem ( 1.1 ) given in the introduction. We assume in this part 



that g does not depend on time. Let a > 0, Afi > and M2 > 0. In the following, we assume that 

^2 : .-j^^j^u,,, — „„j ii„ii _ ^ n.f 



g £ (Fe) is non identically zero and \\g\\ „^ ,^ < Mi, 

H2 (Te) 

q G B'o i{Tq) is such that H^Uss (Fo) ^ M2 and , q > a on Fq. 
Let us prove the following theorem: 



(4.16) 



Theorem 4.16. Let a > 0, Mi > 0, M2 > and uq e H n H^{Qy. We assume that g satisfies (4.151 and 



for j — 1,2, qj satisfies (4.16). We note {uj,pj) the solution of (1.1) associated to qj. Let K be a compact 
subset of {x € Tq Vi ^ 0} where (wi,Ti) is the solution of ( |4.4[ ) with q — qi and let m > be a constant 
such that \vi\ > m on K. Then there exist C{a,m, Mi, M2) > and C'l > such that 



m - q2\\L-^{K) 



< 



In 



C(a,m,Mi,M2) 

Ci 



Ml-M2|lioc(o, + o„;I,2(r^)2,+jtV(«l-n2)-«||j^oo(o, + ^;i2(r,)2)+l|Pl-P2|li^(0, + oo;i,2(r_^,)+||V(pi-p2)-n||_Loo(o^ + ^;_L2(re) 



Remark 4.17. Due the method which relies on semigroup theory, we need measurements during an infinite 
time. 



Proof of Theorem \4.16\ Let j — 1,2 and (vj, tj) be the solution of the stationary problem (2.2) with q =^ qj . 
According to Proposition 2.6 {vj,Tj) belongs to _ff*(f2)^ x H^[Vl) and moreover, thanks to assumptions (4.15 ) 



and (4.16), there exists a constant C{a, Mi, M2) > such that 



+ WtjWhHU) < C{a,Mi,M2). 



(4.17) 



We denote {wj,Trj) = [uj — Vj,pj — tj). Thanks to Theorem 4.3 we are able to estimate \\qi — q2\\L^(K) 
with respect to an increasing function of (ui — t'2)|re ^^"^ ('''1 ~ '''2)|rc ^^'^ their respective gradients in 
norm. Our objective is now to compare the asymptotic behavior of ui — U2 and pi — p2 to the solution of 
the stationary problem vi — V2 and ri — T2. More precisely, we are going to prove that: 

\\wj{t, .)\\m{ny + hjit, ■)\\H^n) < git), 



where g is a function which tends to when t goes to +00. This inequality, combined with Theorem 4.3 
will allow us to conclude the proof of Theorem |4.16[ 

We have that {wj, ttj) is the solution of the following problem: for t > 



dtw — Aw + Vtt 
V • w 
Vw • n — Tm 
Vw • n — Trn + qjw 



0, in n, 

0, in n, 

0, on Fe, 

0, on Fq. 



13 



completed with the initial condition w{0) — uq — Vj. Let t > 0. We have from the theory of analytic 
semigroup that: 



Wj(t,.) 



-tA 



«.«;,(0,.). 



(4.18) 



Let ?7 > 0. There exists a constant C > independent of qj such that: 



(4.19) 



where fj, is given by (4.121 and where || || is the norm operator. We obtain from Proposition 4.13 that: 

\\wj{t,.)\\Hi(n)^ < C{a, M2)\\AgWj{t, .)\\Hi{nr- 



Then, since Wj{t, .) is given by (4.18), and using Corollary 4.11 plus estimates (4.191 with ^7 = f and (4.171, 
it follows: 



\\wjit,. )\\minr < C{a,M2)\\Ale~'^-^w,iO,.)\\mny < Cia, M2)'^i\\uo\\mny + \\vj\\mnr 
< Cia,UQ,Mi,M2)^. 



Using the regularity result for the stationary case given in Proposition 2.6, we have that: 
Note that, thanks to Proposition |4. 13| we have: 

\\dtWj{t,.)\\L2^!^^2 = \\Ag^Wj{t,.)\\L2(^i2-)2. 

Thus, since Wj{t, .) is given by ( 4.18| ), we deduce from estimates (4.19) with 77 = 1 and (4.171 that: 



(4.20) 



(4.21) 



kj(t, .)lk2(n) < C'{a,uo,Mi,M2)- 



-fit 



(4.22) 



Rema rk tha t {uj,p£)€ L'^{0,+oo; H^{n)^) x L°° (0, +(X); H'^{n)). Let ly > 0. In fact, thanks to equa- 
tions ( |4.20[ ) and ( |4.22| ), we obtain that (wj,7rj) e L°°(z^, +00; i/3(17)2) x L°°{iy,+oo; H^{n)) and since 
Uj = Wj + Vj and pj ~ ttj + r,- , we deduce that {uj,pj) G {ly, +00; H'^{i})^) x {i/, +00; H^{i}j). 
Moreover, thanks to Corollary [24 we have iuj,pj) € L°°{0,v;H^{nf) x {Q,v; H'^iVl)). 



We are now able to prove Theorem 4.16 We have from (4.201: 

- W2||L2(r^)2 < C(a, Uq, Ml, Af2)^- + - M2|lL°=(0,+oo;L2(r,)2). 

Then, passing to the limit when t goes to infinity, we get: 

- W2||L2(r^)2 < - M2||L°'(0,+oo;L2(r^)2). 

We prove similarly: 

||V(-Di - V2) ■ n||i2(r^)2 < ||V(wi - U2) ■ ri||L="(o,+oo;L2(re)2)- 
In the same way, but using now (4.221, we obtain: 



||t1 -T2||L2(r^) < Ibl - P2||L°=(0,+oo;L2(rJ), 



and 



||V(ti - T2) ■ n||i2(r^) < \\V{pi -P2) ■ '^IU-(o,+oo;L2(re))- 
To summarize, we have obtained: 

Iki - V2\\l'^(t,Y + ~ ""2) • n\\L-^(T^)2 + ||ti - r2||i2(r^) + |1V(ti - T2) • ri|iL2(rj 

< ll"l-W2|U°=(0,+oo;L2(re)2) + l|V(ui-W2)-"||L~(0,+oo;L2(rJ2) + ||pi-p2||L==(0,+oo;L2(re)) + l|V(Pl-P2)-n||L~ 



Applying Theorem 4.3 to {vj,Tj) for j = 1, 2, we obtain the existence of positive constants C(m, Mi, M2, a) 
and Ci such that 

C(m,Mi,M 2,a) 

Ci 



ki - q2\\L^{K) < 



In 



Ikl— «2|| j,2(r^)2 + ||V('Ui-t)2)-n|| i2(r^)2 + ||Ti-r2|| 1,2 (r^) +11 V(Ti-r2)-n|1^2(p^) 

We conclude by using the fact that the function x h^T) increases on R'^. 



a 
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4.3.3 The flux g depends on t 



We restrict our study to the case where g is cohnear to the outgoing normal n: g — n n. 
Let a > 0, Afi > and M2 > 0. We assume that: 



and 



q e B5 i(ro) is such that WqWbs (To) ^ -^2 and , q > a on To. 

2 ' j ,1 

Let us introduce h such that: 



h ^ H2 (Fe) is non identically zero and H^Hifj (Fq) ^ ^'h 



We suppose that: 



lim ||«(i,.)-/i|L 3 ^ + / e-^(*-«)||at«:(s,.)l 



'H2{r 



H2 (Fe) 



= 0, 



(4.23) 
(4.24) 

(4.25) 
(4.26) 



where fj, is given by equation (4.121 



Theorem 4.18. Let a > 0, Mi > 0, M2 > and uq € H^{n)'^ f] H. We assume that h and k satisfy 



respectively (4.25) and (4.231 and for j = 1,2, Qj satisfies (4.24|. We denote by (uj,pj) the solution of (1.1) 
associated to qj. Let K be a compact subset of {x GTq vi ^ 0} where (ui,ti) is the solution of 

-Au + Vr = 0, infl, 

V • w =0, in D,, 

Vv-n — rn — hn, on Fg, 

\/v-n — Tn + qiv — 0, on Fq. 



We assume that (4.26) is verified. Then there exist C{a,m, Mi, M2) > and Ci > such that 

hi - q2\\L'-(K) 



< 



C{a,m,Mi,M2) 



In 



«l-«2|li«J(o, + oo;t2(r^)2) + l|V(«l-«2)-ri||j^^(0_ + ^.j^2(r^)2) + ||pi-p2|li,^(n, + ocst2(r^))+l|V(pi-p2)-n|Loo(o, + ^;I,2 



Proof of Theorem 4-18 For j — 1, 2, we decompose Uj into Uj = Vj + wj where {vj, Tj) E iJ'*(il)^ x iJ'^(ri) 
is the solution of the stationary problem: 



-Aw + Vr =0, 

V-v =0, 

Vf ■ n ~ Tn = hn, 

Vw • n — TO + qjV — 0, 



and (wj , TTj ) is solution of the following problem: 



dfW — Aw + Vtt = 0, 

V-w =0, 

Vw • n — Tin = (k — h)n, 

■ n ~ im + qjW = 0, 

w{0,x) = uo{x) — Vj{x), 



in Q, 
in Q, 
on Fe, 
on Fq. 



in (0, +00) X fl, 
in (0, +00) X f2, 
on (0, +00) X Fe 
on (0, +00) X Ti, 
in n. 



We would like to perform the same reasoning as in Theorem |4.16[ That is to say, using the fact that we are 
able to estimate \\qi — q2\\L^{K) with respect to an increasing function of (vi — f2)|rc and (ri — T2)|r^ in 
norm, we want to compare the asymptotic behavior of ui — U2 and pi — p2 to the solution of the stationary 
problem vi — V2 and ti — T2. More precisely, we are going to prove that: 



\wj{t, .)\\H3(n)2 + ||V7rj(t, .)||ffi(o)2 < g{t). 
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where 5' is a function which tends to when t goes to +00. Since the function k depends on t, there will be 
one more step than in Theorem 4.16 and that is why we assume (4.26). 

We divide (wj,iTj) into two terms: Wj = vP^ +Wj and TTj = p^j +Pj, where (it^,p^) is solution of 



' dtu° - Au° + \'p° = 0, 



V • 



= 0, 



in (0, +00) X Q, 
in (0, +00) X Q, 



and (wj,pj) is solution of 



Vu" • n — p'^n 


- («- 


h)n, 


on (0, 4 


-00) X Te, 


Vu" • n-p"n + Qju" 


= 0, 




on (0, 4 


-00) X Ti, 


f 


= 0, 




in 0, 




dtw — Aw + Vp = 


0, 




in (0, - 


foo) X fl, 


V -w = 


0, 




in (0, - 


foo) X fl, 


Vw • n — pn = 


0, 




on (0, 


+00) X Te, 


Vw • n — pn + qjW = 


0, 




on (0, 


+00) X Ti, 


iD(0,x) ~ 


uo[x) - 




in n. 





Let t > 0. Using the same arguments as in the previous subsection, we prove that: 



\\wj{t,.)\\H^n)^ < C{a,uo,Mi,M2)- 



-fit 



and 



\\Pjit,.)\\HHn) < C{a,uo,Mi,M2) 



— fit 

IT' 



(4.27) 
(4.28) 



It remains for us to bound ||u°(t, .)||^3(q)2 and ||Vp"(t, .)||//i(o)2. Let t > 0. We are going to prove that 
there exists a constant C{a,Al2) such that: 

\\u°it,.)\\HHn)- + \\^P°it,-)\\HHny 



'H2 (r.) 



ff2 (Fe) 



if 2 (Fe 



ds 



(4.29) 



If inequality (4.29) is satisfied, we can end the proof of Theorem 4.18 



||wi(i, .) - W2{t, .)||//3(o)2 < ||u'j'(t, .) - U2{t, ■)\\m{n)^ + \\wi{t, .) - W2{t, .)\\H3(n)2, 

||V^i(t, .) - V7ri(t, .)\\miny < l|Vp?(t, .) - VpO(t, + l|Vpi(t, .) - \/p2{t, ■)\\minr, 

and in the following two estimates, the right hand side tends to when t goes to infinity thanks to inequal- 
ities ( |4.27p , ( |4.28| ) and assumption ( |4.26| . 

We introduce {yj,Pj) the solution of 

~Ay + Vp = 0, inn, 

V • y =0, infl, 

Vy-n — pn = (k ^ h)n, on Fg, 

\7y ■ 71 — pn + qjy = 0, on F^, 



2.6 



for all t > 0. We know that {yj{t, .),pj{t, .)) e H'^{ny x H'^(n) and satisfies, thanks to Proposition 

■)\\HHny + \\pAt, ■)\\HHn) < C{a, M2)||«:(t, .) - H^.^^^^. (4.30) 

3 

Remark that yj{t, .) belongs to V{Aq.). Indeed, there exists a unique p{t, .) G H'^{Q) solution of 

(4.31) 



Ap 


= 0, 


in ri, 


P 


— K — h, 


on Fe, 


P 


- 0, 


on Ti, 
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for alH > and there exists a constant C > such that 

\\pit,.)\\HHn)<C\\^it,.)-h\\^.^^^^. 

Then [yj , pj + p) is solution of 



(4.32) 



Vj/ • n — (p + p)n 
Vy -n - {p+p)n + qjy 



Vp, in il, 

0, in 17, 

0, on Te, 

0, on To, 



for all t > 0. Remark that, since Vp e L^(ri), we have that yj{t) G 'D{Aq^) by definition of I?(v4g^. ). Notice 
that the fact that g is colinear to n is important here to do the change of variable in the pres sure. We 

1 3 I ' 

deduce from Aq.yj{t) — Vp{t) e V — ^{Aq.) that yj{t) e ^{Aq.). Moreover, using Corollary 
inequality (4.32 1, there exists a constant C > such that: 



4.11 



and 



UlvAt^ ■)\\L^n)- = WAq^yjit, ■)\\m(ny = WVPimminy < CMt, .) - h\\^s^^^^, (4.33) 



that is to say, using moreover (4.30): 



\\y,{t,.)\\^^j^<C{a,M,)Mt,.)-h\\^s 



(4.34) 



We can use the same argument, replacing k ~ h by dtK, to prove that dtyj{t, .) G 'D{Aq.) together 
with the estimate 

\\dtyAt.-)\\^^^l^<C{a,M2)\\dtK{t,.)\\^.^^^^. (4.35) 

Let us consider Wj — vP^ — yj and p^ = p° — pj . The couple {wj , Pj ) is solution of 

9ti« — Aw + Vp = —dtyj, in (0, +oo) X f2, 

V • w =0, in (0, +oo) X n, 

\7wn-pn = 0, on (0, +00) x Fe, (4.36) 

Vw-n — pn + qjW — 0, on (0, +00) x F^, 

w{0,x) = —yj{0,x), in fl. 

We know that Wj is given by: 



Wj{t, .) 



■^yj{0,.)~ e-^'-^^^ydty,{s,.)ds 



Using the family {4'q^)i>i defined by Proposition 



4.10 



have: Wj{t, .) = J2i>i Ci{t)(j)\^, with 



Ci{t)^-e 'S(y^.(O,.),0M^.(^). 



e ^* '^^•'={^tyj[s,.),(|)[^L2(^^Yds. 



Thus, recalling that (A^^.);>i satisfies (4.121, there exists C > such that: 

Ci{tf < 2e-2*''(y,(O,.),0^^)^,(^^, +C /*e-(*-^)'^(a,%(.,.),<,) 



ds. 



We obtain from estimates (4.341 and (4.351: 



Remark that, thanks to Proposition |4.13| and Corollary |4.1H we have: 



(4.37) 



\wj{t,.)\\HHnY < C{a,M2)\\Aw,{t,.)\\Hi(^nr ^ C{a, M2)\\A^Wjit, .)\\L2^nr < C{a,M2)\\wj{t,.)\\ j 

(4.38) 
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To summarize, using (4.38) and (4.37), we obtain the desired estimate: 



Using now the regularity resuh for the stationary problem given in Proposition |2.6[ we have: 

l|Vpj(i, ■)\\HHnr < C{a, Ah) (||9t%(i, ■)\\HHnr- + Hw^it, .)\\HHnr-) ■ 
Since Aq.wj = —dtUj — dtWj, we obtain: 

l!Vp,(t,.)||/fi(n)2 < C(a,M2) {\\dtyj{t,.)\\HHnr + \\Aq,Wj{t, .)\\„i^^)2) . 
Thanks to Corollary 



4.11 



and (4.37), we obtain: 



we know that \\Ag^Wj{t, .)\\Hi{n)2 — \\Aq^Wj{t, .)\\]^2(^fiy2 . Therefore, using (4.35) 



The estimate (4.29) follows from = Wj + yj, = Wj + pj and inequality (4.30). 



□ 



Remark 4.19. Let I e Hf^^{0, +oo; H2 (Fg)) and h e H2 (T^). Assume that there exists 9 > such that: 

supe*^||/(i, .)ll^i(r.) + ll^*^(*'-)ll/fi(r.)) < 



Then k = h + I satisfies (4.26). We note that a particular case of function satisfying (4.26) is given by 



l{t,x) — uj{t)p{x) where uj e Hf^^{0, +00) , p £ H2 (Te) and limt_>oo e uj(t) — limt_>.oo e uj' {t) = 



4.4 Conclusion 



To conclude, we have proved, under some regularity assumptions on the open set SI and on the solution 



{u,p) of system ( 1.1 ), logarithmic stability estimates for the Stokes system with mixed Neumann and Robin 
boundary conditions. Due to the method which relies on the Carleman inequality proved in [3], these 
estimates are valid in dimension 2 and we need measurements on the whole boundary part Fg. 
Our result could be improved in many different ways. In particular, a first concern could be to prove a 
logarithmic stability estimate which is valid in any dimension. It would be also interesting to know if local 
measurements on a part of the boundary included in Fe still allow to get the same result. Moreover, in 
our stability estimates, we need measurements on Fe of u, p, Vit • n and Vp • n, while the identifiability 
result given by Proposition |3.3| only requires information on u and Vu ■ n — pn on Fg. Therefore, it might 
be interesting to know if it is possible to obtain a stability inequality with less measurement terms and in 
particular, if it is possible to get rid of the gradient term Vp • n. 
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A Existence and uniqueness for the unsteady problem 



We study the regularity of solutions of the unsteady problem: 



ut — Au + Vp 


= 0, 


in (0, T) X 0, 


V • u 


= 0, 


in (0, T) X 17, 


Vu • n — pn 


= 9, 


on (0,T) X Te 


Vu ■ n — pn + qu 


= 0, 


on (0, T) X To 




= -"0, 


in r2. 



We are going to prove Theorem |2.7| First of all, as a preliminary result, we prove the following existence 
result: 

Proposition A.l. Let T > 0, M > 0, a > and uq G H . We assume that g e L2(o,T;L2(rg)rf) and that 
q e L°°{Tq) such that q > a on Fq. There exists u G L^{0, T; V) such that for all v € V in the distribution 
sense on (0, T): 

d f f ^ ^ f 

u ■ V + / Vu : Vv + qu ■ V 



dt 



9 ■ V, 



and \fv e V, 



u{0) ■ V — / Uq ■ V 



(A.l) 
(A.2) 



Proof of Proposition A.l . We begin by proving, using a Galerkin method, that there exists u <E L^(0,T; V) 
such that 



Vw e t/,VV' e C^(0,r) such that V^(T) = 



Jn 



u(t, x) ■ v{x)ip' {t)dxdt + 



Jo. 



Vu{t,x) : Vv{x)^{t)dxdt 



q(x)u{t^x) ■ v{x)ip{t)dxdt — I uo{x) ■ v{x)dx 



g{t,x) ■v{x)ip{t)dxdt. (A.3) 



Let (wi)igiM be a Hilbert basis of V which is also an orthogonal basis of H. For each rt e N, we define an 
approximate solution as follows: we search it„ G Vn ~ Vect{wi}i^i^n which satisfies 



/ M„,t -Wj + / Vm„ • Vwj + / qun -Wj = I g ■ Wj^j G {1, . . . , n}, 
Jn Jn Jto Jr^ 

n 

fc=l 

Note that Un,t denotes dtUn- 

Let t e [0,T]. We decompose .) in the Hilbert basis: 

n 



(A.4) 



We denote A= 



Jj-j Wi{x) ■Wj{x)dx 



l<.i,j<.n 



{Ci{t))i<i<n and L{t) — (Jp g{t,x) ■ Wi(x)o?x)i<j<„. We can rewrite the system (A.4) in the form: 

{Aat)+Bm^L{t), 

Since the matrix A is invertible, the system has a unique global solution ^ G H^{Q,TY. We are now going 
to prove that there exists a constant C > independent of n G N such that: 



sup 
o<t<T Jn 



Jn 



|Vu„| 



T 

Jn 



"„r < c. 



(A.5) 
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Let t € [0, T]. Multiplying the first equation of (A. 4) by and summing over j for j — I, ... we obtain: 



"'si Jo Jn Jo JFo Jo Jr^ 



(A.6) 



Let e > 0. Thus: 



„|2 



|2 ^ / / l„|2 I ^ / IL, ||2 



9-Un<C / + / K|^<C/ / \g\' + e h„||^,(^).. 
"'o Jr^ Jo Jr^ Jo Jr^ Jo 

Choosing e small enough and using the fact that q > a on Tq, we obtain: 

sup / iu„p+ r f iv^.„p+ r / \u,,\'<c{r / \g\^+ / kp). 

te[o,T]Jo Jo Jo Jo Jo Jo Jr. Jo 



(A.7) 



This gives (A. 5). According to inequality (A. 5 1, there exists u e L'^{0,T; V) such that, up to a subsequence, 

u„ ^ li in L^{0,T; V). 



Let j G N. Multiplying the first equation of (A.4) by V' e C^([0,T]) such that 'i(j{T) = then integrating 
over (0, T), we get, ^n> j: 

T , .T . 

Un^t{t, x) ■ Wj{x)ip{t)dxdt + / / q{x)un{t, x) ■ Wj{x)tl;{t)dxdt 

Jo JFo 

+ [ [ \/un{t,x):\/wj{x)tlj{t)dxdt=^ [ f g{t, x) ■ Wj{x)^p{t)dxdt. (A. 8) 



Jo 



Jo 



Jr 



Taking into account that: 



Jo 

T 



Un.t{t,x) ■ Wj{x)^p{t)dxdt 



Un{t, x) ■ Wj{x)ip' {t)dxdt — / Un{0,x) ■ Wj{x)tlj{0)dx, 



Jo Jo 



we easily pass to the limit when n goes to infinity in ( A.8). Remark that this inequality is still valid replacing 



Wj by any v & V hy continuity. This ends the proof of the existence of u G L {0, T; V) which satisfies ( A.l ) 
in the distribution sense on (0,T). 



Let us finish the proof of Proposition A.l by proving that the initial condition (A.2) is satisfied 



We deduce from equality (A. 3 1 that 4i{u,v)L2/Q\d G L^(0,T). Consequently, the function t — >• (^(t), f )^2( 



is continuous. This gives a sense to (m(0), w)i2(Q)d. Let ^ G C^(0,T) such that tp{T) = 0. Multiplying (A.l) 
by "0, we obtain: 

{u,v)L2rmd'il;'{t)dt + / aq{u,v)t{j{t)dt = {u{0, .),v)L2my^{0) + / l{v)^{t)dt. 



Comparing with equality (A.3|, we obtain '0(O)(m(O, .) — uq, t')L2(o)<' = 0- Let tp be such that ip{0) ^ 0, we 
have (m(0) — uq, w)i2(Q')d — 0, yv G V . □ 



We are now able to prove Theorem |2.7| 



Proof of Theorem 2.1. We will begin by proving that ut G L^{0,T; H), then we will conclude by using the 
regularity result for the stationary problem from Proposition |2.6[ 



Let t G [0, T]. Multiplying the first equation of (A.4) by and summing over j for j = 1, . . . , n we obtain: 



Wn,tP+ / / qUnUn.t+ / / Vw„ : Vu„,t = / / g-Un,t- 
Jo Jo JTo Jo JO JO JTe 
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We have: 



9 ■ Un.t 



Jt 



Jr, 



9tUn - / 5(0)u„(0)+ / g(t)un{t) 



Let e > 0. Thanks to Cauchy-Schwarz inequaUty and estimate (A. 7), there exists C > 0: 



gtUr. 



^r. 



< C 



^r, 



\9tV+ sup / |5r + hollHi(n)2+e / |V7/„(0|^ 
te[o,T]Jr, 



n 



Using successively integration by parts over (0, T) we finally obtain, choosing e small enough: 



(A.9) 



sup / |V'u„|^ + 
te[o,T] Ja 



"'SI 



\Un,t? < C ||Mo|lHi(o)d + 



"T 



|5tP+ sup / \9\ 



(A.IO) 



te[o,T] JPe 

We deduce that (M„)„eN is bounded ini/i(0,r; iJ)nL°°(0, T; and therefore w G H^{0,T; H)nL°°{0,T] V). 

Then we use th e re gularity result for the stationary problem. For all t e [0, T], we have Ut{t) G 
L'^{flY so by Proposition 2.6 the solution {u{t),p{t)) belongs to H'^{fiy x H^{fl) since the map: 

is linear and continuous. Since {ut,9) e L^{Q,T-L'^{nf) X L2(o,T;i72(re)''), we deduce that (u,p) € 
L\0,T;H'^{n)'^) X L^O,T;H^{n)). 

Let us now prove its uniqueness. Assume that ui and U2 are two solutions and let w = ui — U2- 
Then w e ^1(0, T; 7J) n ^^(o, T; V") and we have for aU veV: 



wt{t) ■v+ Ww{t) : Vw + / qw{t) ■ v ^ 0, w{0) 



= 0. 



(A.ll) 



Taking v = w{t) in (A.ll I, we find: 



1 d 
2dt 



that is to say 



\wit)\^ < / |w(0)p = 0, for aU t e [0,r] 



So ui — U2 on (0, T) X 57. To conclude, thanks to system ( 1.1 1, we obtain pi — P2- 



□ 
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